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Abstract. This paper develops a deeper understanding of the structure and combinato- 
rial significance of the partition function for Hermitian random matrices. The coefficients 
of the large expansion of the logarithm of this partition function, also known as the 
genus expansion, (and its derivatives) are generating functions for a variety of graphical 
enumeration problems. The main results are to prove that these generating functions are 
in fact specific rational functions of a distinguished irrational (algebraic) function, zo(t), 
of the generating function parameter, t. This distinguished function is itself the generat- 
ing function for the Catalan numbers (or generalized Catalan numbers, depending on the 
choice of weight of the parameter). It is also a solution of the inviscid Burgers equation 
for certain initial data. The shock formation, or caustic, of the Burgers characteristic 
solution is directly related to the poles of the rational forms of the generating functions. 

As an intriguing application, one gains new insights into the relation between certain 
derivatives of the genus expansion, in a double-scaling limit, and the asymptotic expansion 
of the first Painleve transcendent. This provides a precise expression of the Painleve as- 
ymptotic coefficients directly in terms of the coefficients of the partial fractions expansion 
of the rational form of the generating functions established in this paper. Moreover, these 
insights point toward a more general program relating the first Painleve hierarchy to the 
higher order structure of the double-scaling limit through the specific rational structure 
of generating functions in the genus expansion. The paper closes with a discussion of the 
relation of this work to recent developments in understanding the asymptotics of graphical 
enumeration. 

As a byproduct, these results also yield new information about the asymptotics of 
recurrence coefficients for orthogonal polynomials with respect to exponential weights, 
the calculation of correlation functions for certain tied random walks on a 1-D lattice, 
and the large time asymptotics of random matrix partition functions. 



1. Introduction 

This paper describes the coefficients in an asymptotic expansion of the logarithm of the 
partition function, Ziy{t), 

log (Zjv(t)/(Zjv(0)) = N^oit) + ei{t) + ^e2(t) + • • • + j^ey{t) + ■■■ 

for certain ensembles of random Hermitian matrices. (See Theorem II. II for a precise state- 
ment.) 

The support of the Natonal Science Foundation, DMS-0808059, is gratefully acknowledged. 
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As is familiar from general contexts of statistics and probability theory, these coefficients 
contain a great deal of information about expectations and correlation functions with re- 
spect to the Gibbs measure that underlies the given random matrix ensemble. In our 
context we are interested in the asymptotic behavior of A^^^ log Zjv(t) (often called the 
free energy) as the size, N, of the matrices becomes large. These ensembles, and hence 
the asymptotic coefficients eg, as well, depend on a coupling parameter t (the analogue of 
the inverse temperature in the classical setting of the canonical ensemble). It is natural 
to view the eg(t) as generating functions for the just-mentioned expectations and correla- 
tion functions and therefore it is of value to know as much as possible about how these 
coefficients depend on the coupling parameters. 

We will derive general closed form expressions for alloi the coefficients, viewed as generating 
functions, in our asymptotic expansions as well as in certain derivatives of these expansions. 
We will show that, with only a couple exceptions, they are all rational functions of a single 
irrational (algebraic) function of t which we denote by zq. More succinctly, we will see 
that, for g > 2, 

eg{t) G Q (zo(— t)) , the function field of zo{—t) over the rational numbers Q. 

It is a standard and fairly elementary fact that a generating function is rational if and 
only if its sequence of coefficients eventually (i.e., beyond a certain order) satisfies a linear 
recurrence relation |20) . The generating functions, such as eg{t), that we will describe are 
not rational in t; however, the result just mentioned suggests that if eg is indeed rational 
in Zq, then it should in fact be possible to recursively construct it from zq. This is, in 
effect, the viewpoint which guided our derivation of the rational form of the eg and other 
generating functions that arise as coefficients in asymptotic expansions closely related to our 
partition function expansions (see Theorems 11.31 and II. 5p . For the fundamental function, 
zo{—t) itself, we are able to describe it very explicitly from several different perspectives. 
For instance, it is itself the generating function for generalized Catalan numbers which are 
familiar constructs in enumerative combinatorics. 

Putting all this together, our results provide a fairly detailed and systematic understanding 
of the asymptotic behavior of correlation functions for the random matrix ensembles that 
we study. While this is interesting in and of itself, there is another level of applications 
for these results. The correlation functions that we are discussing also have a combi- 
natorial interpretation related to diagrammatic expansions that are familiar in statistical 
physics. 

The eg{t) can be topologically identified as generating functions that count labelled maps. 
Roughly speaking, a map is an embedding of an oriented graph into a compact, oriented 
and connected surface X with the requirement that the complement of the graph in the 
surface should be a disjoint union of simply connected open sets. A labelled map is a 
map on which the vertices of the embedded graph have been ordered. (Two such labelled 
maps are equivalent if the respective graphs can be smoothly deformed into one another 
within the surface.) The n^^ Taylor coefficient of eg{t) is an integer equal to the number of 
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distinct equivalence classes of maps to a compact, oriented Riemann surface of genus g that 
can be constructed using exactly n vertices and where each vertex has the same valence. 
(This valence is given by a weight associated to the coupling parameter. See Section 11.21 
and Theorem 11.21 ) The results presented in this paper yield recursion formulae for the 
coefficients of these generating functions, now viewed as counting functions for labelled 
maps, that are based on the partial fractions expansions of the rational forms of the Cg and 
related generating functions. 

Of particular interest is the new light that this application sheds on a 20 year old problem 
stemming from the physics literature on 2D Quantum Gravity [U [El [H] . In mathematical 
terms this problem asks for the detailed description of a particular double- scaling limit 
of the random matrix partition function. We refer the reader to [Hj for the history and 
background. We can give an idea of the contribution of our work to this problem by 
mentioning that ZQ{—t) is the profile of the self-similar solution to an inviscid Burgers 
equation (see Section 12. 3p . From the method of characteristics one may conclude that 
ZQ{—t) has an analytic continuation to the complex t-plane minus a branch cut that extends 
from —oo to a critical value —tc < along the negative real axis. This critical time is in 
fact the shock time at which the Burgers caustic starts to form. One then defines the 
double-scaling limit to be the limit in which t — )• —tc within a sector of the complex plane 
with 3f?t > —tc while at the same time A/^ — )• oo so that the similarity variable 

_ N^/Ht + tc) 
-7 

remains fixed at some finite, non-zero value (7 is a constant to be determined). Because 
of our rationality result, we are, first of all, able to conclude that the coefficients, eg{t), all 
have the same envelope of holomorphy as ZQ^—t), which is the above mentioned slit plane 
with a common singularity at —tc- But, more strikingly, what we find in the double scaling 
limit, is that the dependence on the combinatorial information in ZQ{—t) washes out, so 
that all that remains in the limit is the information that can be deduced from the structure 
of Cg as a rational function of zq, and that rational structure is what we describe in this 
paper. 

One may study the asymptotics of the free energy directly in the double-scaling limit. 
This in fact was done, for a particular ensemble, in \12\ [13j (see also [7]) via a Riemann- 
Hilbert representation of polynomials orthogonal with respect to the exponential weight 
that determines the Gibbs measure for the matrix ensemble under consideration. They find 
that this Riemann-Hilbert problem can be deformed into the Riemann-Hilbert problem for 
the first Painleve transcendent (PI) with respect to the parameter in the double-scaling 
limit. Through this they are, roughly speaking, able to relate the double scaling limit of 
N~'^ log Zjv to an asymptotic expansion of PI in a sector of the plane. 

Our approach is also based on the Riemann-Hilbert problem for orthogonal polynomials, 
but we develop explicit expressions for the double scaling asymptotics of the partition 
function, for a very general class of ensembles, in terms of our representations of the map 
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counting functions as rational functions of zq. This is consistent with the results for the 
special case considered in |12U13j and therefore gives, for the first time as far as we aware, a 
direct and explicit interpretation of the coefficients for the PI asymptotic expansion in terms 
of the graphical enumeration of labelled maps. This connection opens up some interesting 
new directions for investigation that will be briefiy discussed in section 14.51 

The rest of this introduction will amplify upon and detail the themes that have just been 
sketched. 

1.1. Partition Functions for Unitary Ensembles. The ensembles of interest in this 
paper are a sub-family of what are known as the Unitary Ensembles (UE) of random 
matrices; the common space of these ensembles is the space. Tin, of n x n Hermitian 
matrices, M. The Gibbs measures for the ensembles we consider are given by the following 
family of probability measures, 

(1) dfit2^ = exp {-N Tt[V,{M,t2u)]}dM, where 

(2) K(A; t2u) = + t2.A2- 

and the parameter t'2v is assumed to be such that the integral converges. For example, one 
may suppose that 5Rt2i/ > 0. 

The normalization factor Z^\t2u), which serves to make fit be a probability measure, is 
called the partition function of this unitary ensemble. Because the trace, Tr, is invariant 
under conjugation and since M is Hermitian, this partition function reduces to an expres- 
sion which is proportional to the following multivariate density expressed in terms of the 
eigenvalues {Xj} of M [23]. 



(3) 



exp < 



1 " 1 



Xn — Xf 



Z^\t) differs from Z^\t) only by an overall factor which is independent of t. (That is 
because the reduction to eigenvalues comes from a change of variables on M that conjugates 
it to a diagonal matrix; given the invariance of the trace under conjugation, the effect of 
this variable change is only seen in dM and does not involve t.) In all our considerations, 
we will always be working with Z^\t) / Z^j~^\o) . So from this perspective, Z^\t) and 
Z^\t) are equivalent. 

For the coupling weight V^/, we will henceforth often drop the subscript v and just write 
V; similarly, for the coupling parameter t2u we will usually drop the subscript 2z^ and just 
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write t. Which value of i' is intended should be clear from context. (However see section 
11.21 for further comment.) 

In studying asymptotic questions we will always assume that N and n go to infinity together 
such their ratio x = n/N remains constant at a finite, non-zero value. When t = 0, the 
Gibbs measure ([1]) is the measure for the Gaussian Unitary Ensemble (GUE). We will 
typically want to work with a rescaling of our matirx ensembles given by ^ = \/NM and 
which we will refer to as the fine scale. In terms of this calibration, the probability measure 
at t = is 



(4) d^lo{A) = 2-"/27r-"'/2exp{-l/2 TY[A2]}dA 

The general partition function associated to the measure or equivalently the integral 
([3|), can be expressed as an expectation with respect to the GUE ensemble. First observe 
that 



(5) zP{t) = I exp {-N Tr[l/2 +tM^'']}dM 

= iV-"'/2 f exp I - Tr[l/2 + —^^^1 ] dA 

where E^g denotes the expectation with respect to the measure dH). It is then natural to 
define a function that interpolates naturally between the matrix and eigenvalue represen- 
tations of the partition functions: 



(6) rl^it) 



(-p{-]^Tr[A2l}), 



Z7(0) Z7(0) 
where t^^n is called the n*'* scaled tau function. The unsealed tau function is Tn^i. 

A primary motivation for carrying out the detailed asymptotic analysis of the integrals 
([3]) is the insight it yields into current problems of asymptotic combinatorics, specifically 
related to graphical enumeration. The work here builds on the systematic asymptotic 
analysis of ([3]) that was carried out in [9l [inj. In particular we take as our starting point 
the following two theorems from this prior work. 

Theorem 1.1. [9] There is a constant T > such that for (complex) t with ?R.(t) > 
and \t\ < T and for x = n/N in a neighborhood, U, of x = 1, one has a uniformly valid 
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asymptotic expansion 

(7) log7:^jv(*) = n'^eo{x,t) + ei{x,t) + ■\e2{x,t) H . 

as — 7- oo. The meaning of this expansion is: if you keep terms up to order n~'^^ , the 
error term is hounded by Cn~'^^~'^ , where the constant C is independent of x and t in 
the domain {x G U;t £ (0 < Kt) fl {\t\ < T)}. For each £, the function ei{x,t) is a locally 
analytic function of t in a (complex) neighborhood of t = 0. Moreover, the asymptotic 

expansion of t- derivatives of log (^Z^^^ may be calculated via term-by-term differentiation 

of the above series. □ 

Remark 1. In [9], the result was stated only for the special case when x = 1. However, 
it is straightforward to extend this result to the more general setting stated here and this 
was explained in [TO]. For the later parts of this paper (sections 3 and higher) we will take 
X = 1 so that n = N. In this case we will denote the partition function by just Z]\f 
and eg{t) will denote eg{l,t). □ 

Our focus will be on the coefficients, eg{t), which have a direct combinatorial meaning that 
we now briefly explain. The Cg enumerate labeled maps. This has already been described 

in the opening paragraphs of the Introduction. Here we will give the formal and more 
precise definition. 

A map D on a compact, oriented and connected surface X is a pair D = {K{D), [i]) 
where 

• K{D) is a connected 1-complex; 

• [i] is an isotopical class of inclusions i : K{D) — t- X; 

• the complement of K{D) in X is a disjoint union of open cells (faces); 

• the complement of the vertices in K[D) is a disjoint union of open segments (edges). 

We introduce the notion of a g-map which is a map in which the surface X is the closed, 
oriented Riemann surface of genus g and which in addition carries a labeling (ordering) of 
the vertices. 

Theorem 1.2. [9] The coefficients in the asymptotic expansion ^ satisfy the following 
relations. Let g be a nonnegative integer. Then 

(8) eg{t) = Y,\{-ty^g{j) 

in which each of the coefficients Kg{j) is the number of g-maps with j 2v-valent vertices. 

Remark 2. Due to this interpretation, the asymptotic expansion ([7]) is often referred to 
as the genus expansion. 
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Remark 3. It is natural to inquire whether there is a generating function representation 
for unlabeled maps. After ah it is these objects that would seem to have the more purely 
geometric character. Indeed if it were the case that the action of relabeling a labeled map 
always produced a distinct other labeled map, then for each unlabeled map with j vertices, 
one would have j\ distinct labeled maps and '^^jp (the actual Taylor coefficients of eg{t) 
up to a sign) would give the count of unlabeled maps. However, it can happen that a 
relabeling produces the same labeled map and hence this way of counting is not generally 
valid. There are some ad hoc fixes for this but we will not go into these since that would 
take us too far from our main thrust; but, it is worth noting that the kinds of isotropies 
that cause the problem tend to be restricted to low values of the parameters. So for certain 
asymptotic considerations, such as the behavior of these coefficients as j — t- oo, the effect 
of the isotropies "wash out". See Section IT6] for a related discussion. 

One may also ask about the consideration of more general weights, for instance of the form 

^A; ti,...,t„) = iA2+^t^A'", 

m=l 

which would correspond to allowing for odd valences and, moreover, for mixed valences 
when there are two or more non- vanishing coupling constants, tm- In fact versions of both 
Theorems 11.11 and 11.21 were proven in this generality in [9] . The results presented in this 
paper could, in principle, be developed in greater generality along these lines. However, 
to avoid over-complication of the main ideas we will restrict attention here to the case of 
pure (but arbitrary), even valence. □ 

1.2. Different Coupling Weights. There is another (discrete) parameter, z^, in the par- 
tition function ([3]) that determines the exponential weight V{X) in the associated measures. 
(Graphically it determines the uniform vertex valence of the maps enumerated by the eg.) 
The associated coupling parameter t in the weight should more properly be denoted by t2i> 
to distinguish the different weights. However, as is evident from statements in previously 
cited theorems of this introduction, we have in general suppressed the dependence on u in 

(n) 

expressions related to Zj^ or derived from it. In part this is to avoid cumbersome notation. 
In all results and arguments presented in this paper, is a fixed parameter. It does not 
appear in any of the scalings we discuss. The role that u plays should always be clear from 
context. 

But there is a more salient structural reason why it is natural to omit the explicit in- 
dependence: the formulae for coefficients in the genus expansion that we derive all have 
a universal character in u; that is to say, the structure of each of these formulae is that 
of a global meromorphic function of an irrational algebraic function Zo{s) (where s = —t). 
Although V does appear as a discrete parameter in each of these functions, it does not alter 
that function's essential meromorphic form. For example, we will show (see Theorem II. 3 p 
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that, for (7 > 2 




{Z0{S) - l)Qrf(g)(Z0(g)) 



where Q is a polynomial of degree d{g). In particular, the integer- valued functions d(g) 
and o{g) are independent of i^. Moreover, in all cases, the explicit dependence on v in such 
formulae can be completely determined by a finite number of specializations at specific 
values of v {v is a positive integer). Similar structural properties will be shown to hold for 
the coefficients of related asymptotic expansions. 

1.3. The Fundamental Generating Function. The generating functions that we will 
describe are not rational in s; however, as has already been mentioned, we will show that 
they are rational functions of a fundamental generating function zq{s). This breaks the 
problem of describing these generating functions into two parts. The first is to determine 
the rational functions as explicitly as possible to the point that one could say something 
about the associated recursion relation or some other explicit finite description. The main 
point of the present paper is to make progress on this part of the problem. The second 
part is to describe zo{s). Although it is not rational we are actually able, based on [THj, to 
deduce very explicit information about this algebraic function. 

Over the course of this and the next section we will present four different ways of defining 
zo{s), all of which play a role in describing and applying the structure of the coefficients 
of the genus expansion and related geometric expansions. These four different characteri- 
zations are 

(i) as a generating function for the generalized Catalan numbers (section 11.4. ip ; 

(ii) as the leading order term in the large expansion of the recurrence coefficients for 
monic orthogonal polynomials with respect to the exponential weight exp{—NV{X)) 
(section 11.4.4]) : 

(iii) as the self-similar solution of an inviscid Burgers equation (section 12. 3p ; 

(iv) in terms of the endpoints of the support for the equilibrium measure of the large 

limit of random Hermitian matrices (section 12. Sp . 

Based on the previous paragraphs, there are several ways in which we can talk about the 
generating functions in the genus expansion. When we want to think of them directly as 
functions of the coupling parameter we will write eg{t) or equivalently eg{—s) (assuming, 
when we do this, that x = 1). When thinking of its structure as a function of zq we will write 
eg(-2o) or, sometimes, eg{z). Finally the two perspectives are related by eg{—s) = eg{zo{s)). 
Similar conventions will apply to related expansion coefficients such as the Zg first described 
in section 11.4.41 
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1.4. Statement of Results. The first steps toward elucidating the structure of the gen- 
erating functions eg{t) were taken in |10j. There, differential equations for these functions 
were derived that made it in principle possible to recursively solve for the eg in terms of a 
distinguished combinatorial generating function denoted zq{s). The differential equations 
in question will be reviewed in section [2l but we want first to give a brief overview of the 
results of this paper and how they depend on prior work. 



1.4.1. Generalized Catalan Numbers. We begin with our first description of the fundamen- 
tal analytic function zq{s). 



(9) zo{s) = ^ciCjS^ where 

j>o 

c. = 2.(';_-^l)=(.+ l)(^2.^^ and 

^ J ( J-1 ) " + i ( / )■ 

(Section 12.31 explains where this expansion comes from.) When u = 2, (j is the j'*'' Catalan 
number. For general these are the generalized Catalan numbers which play a role in a 
wide variety of enumerative combinatorial problems. For a discussion of these applications 
see [21]. 



1.4.2. Background for the First Main Result. In [T^ the first few coefficients of d?]) were 
calculated explicitly: 

eo{zo) = vizo - 'i-){zo -r)+'^ iog{zo) where 

and 



4i/(z^ + 1) 
3(z^ + l) 
u-1 ' 



61(2:0) = -^log{i/ - {u - l)zo) 



{ZQ - l)Q(zo) 
[i^ - {i^ - l)zof 



62(^0) = ; — where 



^(■^o) = {{-1^' + + 8u'') + {-u^ + 411.3 _ 24^4 _ -^g^5)^^ 

zooU 

+ {AAv - + Mv'^ - llv^ + 8v^)zl + (-12 - 12z. + lOSi/^ - 1321/^ + ASv^)zl 
+ (-12 + 48i/-72z^2 + 48z^3-12z^*^)z^}. 
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Remark 4. The expression for 62 (-Zq) given in [lOj contained additional terms involving 
possible constants of integration. The authors were able to rule out the presence of these 
other terms for most, but not all, values of v. However, it is a consequence of the results 
presented here (Theorem II .Sp that these terms are absent for all values of u and hence the 
expression for 62(^0) given above is universally valid. □ 

There are several striking features common to these first few coefficients. 

First, they are all abelian functions of the single variable zq. (An abelian function is a 
function that can be represented as an iterated path integral of a rational function. The 
rational functions are a minimal subclass of the class of abelian functions.) 

Second, the singularities of these functions are restricted to zq = 0, i^/ii' — 1). 

Third, the coefficients of these generating functions can be determined by substitution of 
the generating function for the higher Catalan numbers, zo{s), into these abelian func- 
tions. 

Finally, although the coefficients of these representations do depend (rationally over Q) on 
u, they have a universal character; in particular, one should note that knowledge of the 
coefficients for relatively few values of suffices to completely determine these functions 
for all values of u. 

The four above-stated features were established in [10] (see Theorem 12.41 of this paper) to 
hold for all the genus expansion coefficients. But we will show that the coefficients, eg, are 
even more tightly structured than these properties might suggest: 

1.4.3. Statement of the First Main Result. 

Theorem 1.3. For g > 2 the coefficients of ^ are rational functions of zq of the form 

[zq - l)Qd{g){zo) 

where Q is a polynomial of degree d{g) whose coefficients are rational functions of u over 
the rational numbers Q. The exponent o[g) and the degree d{g) are non-negative integers 
to he determined. 

Once the rational expression for Cg has been determined, the map counts Kg{n) may be 
read off directly by methods of classical function theory. This will be demonstrated for 
a related counting problem in Section 14.11 More could undoubtedly be said about the 
structure of the Cg, but that will be developed elsewhere. 



eJzo) 
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1.4.4. Recurrence Coefficients. We turn next to a related map enumeration problem which 
also has applications to the theory of orthogonal polynomials with exponential weights. 
Let '7r„^Ar(A) denote the n*'^ monic orthogonal polynomial with respect to the exponential 
weight exp(— A^y(A)). These polynomials satisfy a three-term recurrence relation of the 
form 

Further background on these orthogonal polynomials and their relation to the partition 
functions ([3]) is presented in the Appendix. 

In [10] it was shown that, as a consequence of Theorem 1 1.1^ 

Theorem 1.4. [lOj The recurrence coefficients for the monic orthogonal polynomials with 
weight exp(— A^y(A)) have a full asymptotic expansion, uniformly valid for t € (0 < 3?t) H 
< T) and x = n/N in a neighborhood, U , of x = 1, of the form 

(10) bl^{t) = x{zo{s) + ^zi{s) + ^Z2{s) + ---) 

where s = —x'^~^t (see Section \2.1.1\ for the explanation of this scaling) and 

(11) bl^{Q) = n/N = x. 

Remark 5. This result gives us our second interpretation of zq: it is the leading order 
term in the large N expansion of the recurrence coefficient 6^^(t). The fact that this 
leading coefficient does indeed coincide with ([9]) follows from the derivation that will be 
reviewed in Section 12.31 □ 



The relation of these expansion coefficients to map enumeration is given by the following 

m- 

d-^ z 

(12) -Zo"'^(0) = ^ |s=o =* {two-legged (7-maps with j 2z/-valent vertices }. 

^ ds^ 

(A leg is an edge emerging from a univalent vertex; so that the leg is the only edge incident 
to that vertex. Thus the maps being counted here have exactly two vertices of valence one; 
all other vertices have valence 2v.) 
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1.4.5. Background for the Second Main Result. In [lOj . the first few higher order terms in 
the expansion (|10p were calculated exphcitly: 



(13) zi{zo) = 4 J —-4 ^ 

r i/(zy + 2)/12 -u{^u + 2)/12 1/2/6 
= ^0 S : — + : — 3 + 



(i/ - (i/ - l)zo) (z/ - (z/ - 1)2:0) (i^ - - 1)^0) 

Zojzo - l)-P4(^o ) 

1/ - (z/ - 1)20) 



(14) 2:2(2:0) = — — TTT^ where 



^4(^0) = Y^(^ - 1)^ [(2^' - 14^^' + 241^') 
+ (-12z/^ + 148Z/'' - 5A6u^ + 758z/^ - 252z/'^ - 96i/^)zo 
+ (264z/2 - I510u^ + 25551z/^ - 500u^ - 1789z/^ + 8407/"^ + lUv^)z'^ 
+ (-5361/ + 1396Z/2 + 912i/3 - 4596z/'' + 2492z/^ + 1296i/*^ - 868z/^ - 96z/^)zg 
+ (168 + 234z/ - 1467i/2 + 558u^ + 19022/"^ - 1446i/^ - 267u^ + 2941/^^ + 24u^)z^] 

(15 23(^0 = 7 7 .X where 

Pj^zq) is a polynomial of degree 7 (for the explicit expression, see pp. 66-67). 



1.4.6. Statement of the Second Main Result. We will show that the pattern which appears 
to be emerging in the preceding examples is in fact the general structure for the coefficients 
Zg. Note that we adopt the same notational conventions for functional dependence of the 
Zg that were described for the Cg at the end of section 11.31 

Theorem 1.5. 

/ s _ 2:0 (zq - 1)^3^-2(2:0) 
" {u-{iy-l)zo)^9-l 



^^^^ - (z/ - l)zo)2» ^ (1/ - (Z/ - 1)20)29+1 + ■ ■ ■ + (z. - (i. - l)zof9-^ j ' 

where Pzg-2 is a polynomial of degree 3g — 2 in zq whose coefficients are rational functions 
of v over the rational numbers Q. 



This theorem has a number of remarkable corollaries. 
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1.4.7. Residue Formulae for Graphical Enumeration. For instance, we find that 

* {two-legged g-maps with j 2i/-valent vertices } 
3p-i+i 

= ci ^ a^^'^\v) where 

1=0 

= [{j-l)y-{2g + l + {j-2))]a^t'''\^) + ^['^9 + ^+i3-'i)]a'^^X'\^) and 
a^^'^\v) = a^^\v), which is imphcitly defined as the residue in (jl6p . 

1.4.8. Relation with the First Painleve Transcendent. When = 2 we will show that 

4?_i = -2^^-1(2/3)^/2^, for 5>1 

(0) X 

a} = Oifl, 

where the {a,} are the coefficients in the asymptotic expansion of a class of solutions, ()70p . 
to the first Painleve equation. As will be discussed in section 14. 5^ this points to a broader 
connection of these enumerative coefficients, for general u, to the first Painleve hierarchy 
and of the other fundamental residues, a^f\i'), in (flB]) . to the higher order terms in the 
expansion of the double scaling limit of the recurrence coefficients. 

1.5. Outline. In Section [2] we will review the differential equations for the generating 
functions Zg{s) and also provide additional characterizations of the fundamental generating 
function zq(s). In particular we will describe the relation between the polar singularity of 
Zg{zo), at zq = vjiy — 1), and caustic formation in a Burgers equation for zq. With this 
background, we give the proof of Theorem 11.51 in Section [3l The differential equations for 
the generating functions Cg are reviewed in Section 14.21 and then Theorem 11.31 is proved. 
Applications of and corollaries to Theorem 11.51 are presented in Section HI These include 
those results mentioned at the end of the previous subsection. Finally, in Section 14.61 we 
draw some analogies with and possible connections to other recent work in the theory of 
map enumeration. 

Appendix A presents the essential elements that we use related to orthogonal polynomials 
with exponential weights, their recurrence coefficients and the Toda Lattice equations which 
describe how these recurrence coefficients transform when the parameters in the weights 
are varied. We also indicate how the differential equations reviewed in section [2] are related 
to a continuum limit of the Toda lattice. In Section IA.3I a new combinatorial formula is 
derived for the forcing coefficients that appear in the differential equations for Zg. These 
have relevance for calculating the correlation functions for certain tied random walks on a 
one-dimensional lattice. 

In Appendix B we briefly review the one-point correlation functions for eigenvalues of the 
UE ensembles and the Riemann-Hilbert methods that were used to deduce the structure 
of their asymptotic expansions. The purpose of this is to give points of reference for 
some of the theorems quoted in the introduction (particularly Theorems 11.11 and II. 4p . 
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The Appendix also outlines an extension of these prior results to describe the large time 
behavior of the coefficients Zg and eg. This extension is needed in the proofs of Theorems 
[L5]and[L3l 

2. Background on Differential Equations for Higher Genus Coefficients 

A complete differential characterization of the asymptotic coefficients, Zg, appearing in (jlOp 
is given by the first three theorems below. This characterization is founded on a continuum 
limit of the well-known system of differential equations (the Toda Lattice equations) which, 
in this setting, describe how the recurrence coefficients evolve with t. In Appendix A the 
reader will find the pertinent background on Toda equations and a summary of how their 
continuum limit (j22p is derived. 

2.1. Scalings and Notation. The principal results of this paper all involve studying 
the behavior of certain functions of the partition function in the asymptotic limit as the 
parameters N,n and/or t tend to infinity (or some critical value), typically with some 
combination of these parameters being held fixed. In this subsection we summarize the 
various scalings that will be considered. We will also take this opportunity to state some 
notational conventions that will be used systematically throughout the paper. 

2.1.1. Fine Scaling. For several reasons, and in particular when we discuss, in section [2] the 
derivation of the differential equations which yield explicit expressions for eg and related 
generating functions, it will be most effective to use the fine scaling representation of the 
partition function in terms of the GUE expectation that was given in ([6]) . There the natural 
scaled variable to consider was jp^- There is a closely related fine scaling that will play 
a role in our second main result related to recurrence coefficients (see section ll.4.4p . The 
scaling variable in this context will be given by jLi , where s serves as a similarity variable 
relating x and t (see (jlSp ). Recall that n and grow at the same rate so that x = n/N 
remains fixed at a finite, non-zero value. The relations between all these scalings may be 
succinctly summarized as follows: 



(18) s = -x^-H. 

The macroscopic variable will serve to relate the partition function to standard expres- 
sions for tau functions and Hirota formulae that play a role in the developments described 
in sections 11.4.41 and [2l Note that when x = l,s = —t. The minus sign between s and t 
has been introduced because, although the current usage of t is standard for most repre- 
sentations of random matrix partition functions, s (= —t) is the natural variable for the 
generating functions that appear in the genus expansion. See, for instance, formula ([8]) 
above. The factors of ^ are introduced in (jl7p in order to remove an overall factor of ^ in 
the continuum Toda equations (j22p . 
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2.1.2. Spatial Fine Scaling. We will in addition sometimes need to compare asymptotic 
expansions in {n"*"} to those in {{n + i)^'^'} where i is small (and bounded) in comparison 
to n. Of course, these are equivalent asymptotic gauges and it is natural to relate them by 
a power of the spatial scaling variable 

(19) w = (l + - 

y n 

which we think of as a quantity close to 1. This is similar to x but arises in a different 
context. This variable also serves to mediate between the associated fine scaling variables. 
For instance if 

(20) = -'^0= ' then, 

(21) ~s = w'^-^s. 

2.2. Continuum Limits and Cluster Expansions. 

Theorem 2.1. [lOj The continuum limit of the Toda Lattice equations l^94\ ) with respect 
to the time scaling jj?] ) and the lattice scaling il9\) as n ^ oo is given by the following 
infinite order partial differential equation for f{s,w): 

/, = fM (n-i;/,/^,...,/^^,...) = 

(22) Ci/f fw H 9-^1 (/) fwj fwwi fwww ) + ■■■ 

H fw, /tt,(2), • • • , /u,(29+l)) H 

for {s,w) near (0, 1) and initial data given by f{0,w) = w. 

(23) i^f= E nSlv^^^'^'^^'nf^) 

X:\X\=2g+lBl(X)<u+1^^3'^^^^>- j ^ ^' ^ 

where A = (Ai,A2,...) is a partition, with Ai > A2 > A3 > of 2g + 1; rj{X) = 
#{Ai|Aj = j}; £{X) = ^■^ length of X; and \X\ = Yli^i ^i^^ of X; and 

d^^'^^ are coefficients to be described in the next proposition. 

The following proposition gives an explicit closed form expression for the coefficients d^^'^\ 
This is a new result. Its proof will be given in lA.3[ 

Proposition 2.2. 

d'^x^^ = E 2 mA (/Ui - 7?i, . . . , jiy+i - r]y+i) 

(1/ + 1, zv, . . . , 2, 1) C /X C (2i/, 2u-l,...,v) 
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where TZ is the set of restricted partitions (meaning that /ii > ^2 > • • • > fJ-u+i), {vit ■ ■ ^ "Hv+i] 
{21/, 2i/ — 2, . . . , 2, 0), and mx{xi, . . . , x^+i) is the monomial symmetric polynomial associ- 
ated to A [2Tj. The relation of inclusion between partitions, p ^ means that > pj for 
ah j. 

Remark 6. As a polynomial in the ^-derivatives of /, Fg'^^ is reminiscent of the partial Bell 
polynomials, B^x\,e{X)'i however, it is different from these in that it depends on an additional 

parameter, v + 1, and the coefficients, which depend on d^^'^\ are combinatorially more 
complex. We note also that these coefficients differ slightly from those presented in |10] 
in that the analogous coefficients in the former paper incorporated the parts multiplicities 
factor, l/Hj^i') i^^^o coefficient itself. 

In common parlance, one would take the continuum limit in Theorem 12. II to be the limit of 
these equations when n — )• 00; i.e., the leading order terms which constitute a finite order 
(in fact first order) pde. We will in fact be focussing on this leading order equation in the 
next two subsections. However, for the general analysis we want to carry out in this paper, 
we will need to look at the successive higher order corrections present in (I22p . In fact we 
will need to consider these corrections to all orders. So we are interested in the hierarchy 
of finite order pdes generated by this expansion. This is completely in the spirit of the 
asymptotic analysis of weakly nonlinear differential equations. As in that setting, we posit 
the solution to have the form of an asymptotic expansion: 



(24) 



fis,w) 



1 



1 



fo{s,w) + z^lf^^^^'^) ^ Zl^-^s^^'"^) + 



n 



n 



The hierarchy of equations in question may then be developed by inserting (|24|) into 
and then successively, as g increases, collecting terms of order in n~^^. It should be clear 
that at each order these equations will be partial differential equations, but in several depen- 
dent variables: /o, /i, . . . , fg-i, fg- It should also be clear that this hierarchy is a triangular 
system in these dependent variables: at order 2g the pde is an equation for fg{s, w) forced 
by differential terms in /o, fi, . . . , fg-i which should already be known functions having 
been determined by solving the lower order equations in the hierarchy. These statements 
and the structure of the hierarchy are made precise in the next theorem. 

Theorem 2.3. [lOi The order 2g equation, for g > 0, in the asymptotic expansion of 



on 



(25) 



(26) 



IS 



^ = {{foYifgU + ^{kr-\h)u,fg) + Forcing^, where 



( 



Forcing^ = 



+ 



< kj < g 
\ ki-\ h k^+i = g 



Ft\2g-2] + F^''\2g-A] 
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F^'^\2r] denotes the coefficient ofn~^^ in F^^\ Note, in particular, that such a term cannot 
involve fk if k > r . In ^26\) it is natural, and it will be convenient, to denote the first group 
of terms by FQ'^\2g]. Also, for g ^ 0,fg{0,w) = and, at leading order, /o(0,ri;) = w. 

The structure of (f25]) and (1261) suggests that these equations may be solved recursively 
unless a resonance arises at some stage that obstructs the solution of the linear equation 
(I25p and requires the introduction of a solvability condition in order to continue. Absent 
such resonances, one expects to find expressions for each fg, for g > 0, as functions, in 
some appropriate function class, of just /q. Generally, in the asymptotic analysis of weakly 
nonlinear differential equations, one does not expect to get away without running into 
resonances. However, in our special case we do in fact avoid resonances at all stages. This 
was established in [10]. In fact more was shown there: the formal series (j24p that can 
be built from these recursively constructed coefficients is in fact an asymptotic expansion 
uniformly valid for s negative and near 0. The precise statements are given in the next 
theorem. Appendix IA.2I outlines the essential elements underlying these results. 

Theorem 2.4. [10] 

fg{s,w) = w'-^3Zg{w''-h) 

1 1 

/(S,l) = Zo{s) + ^Zl{s) + ^Z2{s) -\ . 

(ii) The coefficient Zg is an abelian function of zq with singularities only possible at 
zq = and zq = ^/{v — 1). 

(iii) The coefficient Zg is more explicitly presented as a function of zq through the fol- 
lowing integral equation: 

t , ( ( zojs)'^'-'^ f (t^-(z.-l)z) ^ . , . , 

= zg{z,{s)) = ,^,.+3-2, Forcmg,(z)&, 

where 

Forcingg(zo) = Forcingg|^^_^ . 

The terms in Forcingg|^„=i depend only on Zk{zQ),k < g and their derivatives and 
these derivatives can in turn be re-expressed in terms of Zk[zQ) for k < g. (See 
sections \3.1.^ and \3.1.^) . 

Note that (see pO|) ) and xf{s, 1) possess the same asymptotic expansion. Note also 
that since f{0,u>) = w, and consistent with (fTT]l . 

zo(0) = 1 

Zg{0) = for g > 1. 
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A principal result of this paper is to refine part (ii) of the previous theorem to state that Zg 
is a rational function of zq for g > 1. The rational functions are a subclass of the abelian 
functions. Another way to say this is that we establish solvability, without resonances, of 
the differential hierarchy (j25p . beyond g = 1, within the function class of rational functions. 
There is a solvability condition at g = 1 which is to choose a branch of the logarithm at 

Zq = u/{v - 1). 

2.3. Leading Order. The leading equation in the PDE scheme ()22p is 

(28) fs = c,ru 

f{0,w) = w, 

which is an instance of the inviscid Burgers equation, well-known from the theory of shock 
waves. At leading order in n this becomes 

^ = Cu{foy{fo)w 

Substituting the self-similar scaling (p7|) and then evaluating at = 1 it further reduces 
to the ODE 

(29) 4(s) = c^zoisr {zois) + (i/ - 1)sz'q{s)) 
with initial condition zq{0) = 1. The solution is given implicitly by 

(30) I = Zois) - Cuszois)" , 

which is our third interpretation of zq and the one from which our first one, ([9]), was 
derived. 

Remark 7. We mention that, in the case u = 2, (j30p has a relation to the following 
discrete equation that played a key role in the early studies of 2D quantum gravity in the 
physics literature. 

(31) ^KNit) (d,7v(i) + blNit) + bl+iM*)) + blAt) = -. 

In the literature on orthogonal polynomials, this relation is one of the basic examples of 
what are known as Freud equations [22]. In the physics literature it was referred to as the 
string equation. It is an example of a discrete Painleve equation and, in an appropriate 
sense, is a discretization of the ordinary differential equation that defines the (continuous) 
first Painleve transcendent [12^ I13j. With respect to the continuum limit described in the 
previous subsection and based on Section [A^ ([3T1) limits at leading order to ([30|) evaluated 
at v = 2. 

2.4. The Burgers Caustic. As illustrated in Figure [U the solution, zq develops a caustic 
where z'q becomes undefined. (This infinite slope corresponds to the point, in unsealed 
variables, where characteristics for (I28p would first focus or shock.) 
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Figure 1. Burgers Caustic 



From (I29p . and using (130p . one can express the derivative purely in terms of zq as 



(32) z'o = - 

From this one sees that the shock coincides with the proposed unique pole location, 

v -{v - l)zo = 0, 

in the rational expressions for the higher generating functions. For future use we want to 
calculate here the shock time Sc which corresponds to zq = ^/{ly — 1); 



s zo 



u- 1 

zq-I, 



CitZ, 



V 1^0= — 



(33) = ^ ' 



From ()32p it is straightforward to derive the leading local behavior near the caustic, 



(34) ("°-^) = -y7^^(^-^^){i+o(^-^^)}- 

In the special case oiv = 2^ Sc = and one can solve globally for zq from ([30|) using the 
quadratic formula, 

1 - ^/Y^^Ws 
= 24^ • 

Remark 8. It follows from these observations and Theorems 11.31 and ll.Sl that the maximal 
domain of holomorphy, containing t = 0, for both eg and Zg is the set of t S C\{—Sc, oo). 
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2.5. Another Interpretation of zq. The fourth and final interpretation of zq is in terms 
of the spectral density for random Hermitean matrices with respect to the probability 
distribution ([1]). For the class of weights, V^, we are considering and for all t > (see 
section lB.ip . the mean density of eigenvalues for this ensemble limits to a continuous density 
supported on a finite interval [—f3,f3]. This density defines a probability measure, usually 
referred to as the equilibrium measure, which has the following characterization. 

Theorem 2.5. [5] The equilibrium measure fiy is absolutely continuous with respect to 
Lebesgue measure, and 

(35) V(A) = ^X(-/3,/3)(A)v/(A + /3)(/3-A)MA), 

where h{X) is a polynomial of degree 2v — 2, which is strictly positive on the interval [—(3, /3] . 
The endpoints of support ±/3(t) vary smoothly with t (or equivalently s) and 

(36) zois) = ^tf). 

The fact that this characterization of zq coincides with the other definitions given in this 
and the previous section, follows from the requirement that the equilibrium measure should 
be a probability measure. This is an algebraic condition on the endpoints of support that 
is equivalent to (|30p . A detailed explanation of this relation is given in |1U) . 

The equilibrium measure is the unique solution of the following variational problem for 
logarithmic potentials in the presence of an external field [25]. (In this setting, V is the 
external field.) 



(37) sup <^ - / V{X)dn{X) + / log |A - T]\dn{X)dfi{ri) 

where A is the set of all positive Borel measures on the real axis with unit mass. This 
variational problem, and hence the definition of ^uy, can be extended to t G (— Sc,0) [7]. It 
follows that this interpretation of zq extends as well to this negative interval and coincides, 
there, with its other interpretations given in this section. 



3. Rationality of the Asymptotic Recurrence Coefficients 



The first step toward proving Theorem 11.51 is to show that Zg must be a rational function 
of Zq . Given parts (ii) and (iii) of Theorem 12. 4| it suffices to show that the integrand in 
Theorem 12.41 (iii) has no residues at z = and z = vjiy — 1). In |10) this vanishing of 
residues was established for z\,Z2 and 23 (in fact. Theorem 11.51 was fully established in 
these three cases by direct calculation). If there is a term of Forcing^ which generates 
a non-zero residue in the integrand of (2) we call that term resonant. We will show by 
induction that for all g, Forcing^ can be decomposed as a sum of terms each of which is 
non-resonant. 
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Proposition 3.1. (i) Zg is regular at zq = and in fact vanishes at least simply there. 

(ii) Forcing^ may be written as a sum of terms each of which has the form 

/o ^ (/fci)ui(Ji) (•^fc2)ui(J2) ■ ■ ■ (/A;m)ui(j™) 
where < m < + 1 and Ji + • • • + Jm = 2(5 — ki — k2 — ■ ■ ■ — km) + 1, with ki < g. 

(iii) For {k,n) 7^ (0,0), and assuming that Theorem 11.51 holds for all < A; < 5, 

r Ak,n) Ak,n) >| 



(1/ - (Z^ - l)zo)2'=+" (1/ - (z. - l)zo)5^+2'^-l 



The proof of the rational form (jl6p for follows fairly straightforwardly from this propo- 
sition. 

Proof (of Theorem ll.Sp . Inserting expansions of the form (iii), from Proposition 13. 11 for 
each factor of the product appearing in (ii), from Proposition 13. 11 one derives the rational 
form (as a function of zq) of each term in Forcing^ (see Section 13.1.21 below for further 
explanation). One may then deduce, from Theorem 12 .4( 111) . that Zg is a sum of terms, with 
K = kl + k2^ h fern and J = ji + J2 H 1- im, of the form 

V' r (z. - (z. - l)z) f Do ^ 



l)zo J I CyZ''+^-'^9 _ (i, _ l)2)2ff+l 

—m—1 



{V-(V- \)zfK+2J-ra 



dz 



^0 / _^29-2 I ^0 ^ 



D3K+J 

+ (i. - (Z. - l)z)5K+2J-m-l j 



c,{u-{u-l)zo) Ji V J \{u-{iy-l)zY9 

D3K+J —m—1 1 , 



^0 



f -E-O Er^K+2J-m-i 1 , 



C^(zy - (1/ - l)zo) 1 (1/ - (1/ - l)z)2 ' {^u-{u- l)zfK+2.J-m-l 



1 I -E'-l , Eq E^k+2J- 



(38)= ... .^■,^ + ... ... ^2 + 



m— 3 



23-2 1 _ _ 1)^^) {^u-{v- l)zo)2 (1/ - (1/ - l)zo)5^+2J-m-l 



where the leading pole term, with coefficient -E-i, arises from the lower bound of integration 
which insures that this expression vanishes at = 1. The integrand in the penultimate 
line above is clearly non-resonant and it therefore follows inductively that Zg is rational as 
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a function of zq. However, in order to establish the form for Zg stated in equation (jl6p and 
to complete the induction we need to make further use of this result. The total expression 
for Zg is comprised of a sum of series of the form (|38p . each series coming from one of the 
terms into which Forcing^ may be decomposed (see Proposition 13 . 1 ( ii) ) . 

We determine next the maximal pole order of this total expression. The pole order of each 
component series as given by the formula above, is 5K + 2J — m — 1. By Proposition 13.11 
(ii) , there is one constraint, J = 2{g — K) + 1, and two boundary conditions, < m < + 1 
and ^ < K < g. It is straightforward to check that this linear problem is maximized over 
integer values by (K,m) = {g,2) or {K,m) = {g — 1,1), so that, respectively, J = 1 or 
J = 3. Hence the maximal pole order can be at most 5g — 1. 

When J = 1, the terms that contribute to the maximal pole are of the form 



9-1 

(39) l^cJl^-'Y.fmifg-ml^ 

m=l 

as can be seen from (j22p . 

When J = 3, a term yielding this maximal pole order could only come from theFi"^ [2g-2] 
component of Forcing^ as can be seen from ()26p . Moreover, since m = 1, the only term 
within this component giving this maximal order comes from the summand associated to 
the partition A = (3). 

More will be said later about the structure of the coefficient of the highest order pole 
(Section 14. 4p . For now it will suffice to conclude that Zg has the form 

(40) .. = 4r.l7-^^^TrT + 7-^:^^Tr^ + -+ ""''-^ 



{u-{iy-l)zo) {u-{u-l)zo)^ {iy-{u-l)zof9-i 



•^9 — 2g-2 
^0 

for all iy>2. 

By Proposition 13.1 ( i) . Zg vanishes at zq = 0. For this to be possible, the Laurent series in 
0|) must have a factor of Zq^ ^ and hence this expression for Zg becomes 

f n^a) Aa) J9) ^ 

^0 < 7- 7- 77-^977 + 7- 7- n , n9.4-i ^ 



(41) 



Zo{zo - l)P3g^2izo) 



where the factor of — 1 in the last line is inferred from the fact that each term contributing 
to the final expression vanishes at zq = 1. This completes the induction and the proof of 
Theorem 11.51 □ 



3.1. Proof of Proposition [STTl 
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3.1.1. Regularity at zq = 0. Recall from ([32]) that zo{s) satisfies the first order ode 

It follows from this dynamical representation that if zq is initialized at 1, then its trajectory, 
as s — 7- — oo, limits monotonically and asymptotically to the fixed point at zq = 0. Hence, 
to study the behavior of Zg as zq — ^ 0, it suffices to consider the limit of Zg{s) as s — >• —oo. 
By Theorem lB.3l one has an asymptotic expansion for the orthogonal polynomial recurrence 
coefficients that is uniformly valid for s in a neighborhood of the half-line (— c«,0) in the 
complex s-plane, so that 

(42) \bl^^{^-s)-z^{s)-N~^z,[s) iV-'%(5)| < 

where C is independent of s in a neighborhood of the negative s-axis. We proceed induc- 
tively and so can assume that Zj for j < g vanishes as s — )■ — oo (we know, to begin with, 
that zq vanishes in this limit). If one restricts to s on the negative real half-line, then (j42p 
may be re-expressed as 

(43) < N'^b%^j^{-s)-zo{s)-N~'zi{s) N-^'^+^Zg^,{s)) and 

(44) zg + -^ > N'^{b%^^{-s)-zo{s)-N-'z,{s) N-'^+'zg.,{s)) . 



iV2 

We now choose and fix N so that < f • Next we note that 



b%M-^) < ^/^A^P^(^,A)e-^^(^-^)dA 
— )• as s — )• — oo, 

where Pn is the n*^ orthonormal polynomial with respect to the exponential weight exp {—NV{s. 
exp (iV(iA2 - sA^^)). The inequality is a direct consequence of the recursion formula for 
these orthonormal polynomials. The limit in the second line follows because the family of 
densities {p%{s, A)e~^^(''''*')} constitutes a delta sequence as s — >• — oo. 

As a consequence of this observation and the induction, one may choose |s| sufficiently 
large so that the right-hand side of ()43p is less than e/2 and the right-hand side of ()44p is 
greater than —e/2. Therefore, 

for s sufficiently negative. Since e was arbitrary, this establishes statement (i) of Proposition 
\3A[ 



3.1.2. The Structure of Forcing g. By (p6|) . Forcing^ is a graded sum 

Ft'>[2g]+Fi''\2g-2] + ... + Fg[0] 
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where [2(7 — 2r] is the coefficient of n~'^^~^'^''' in . The terms in Fr are in 1 — 1 
correspondence with the partitions of 2r + 1 as described in (I23p : i.e., the terms in this 
component are each proportional to an expression of the form 

(45) fo^™'^^(.fki)w{n){fk2)w(n) ■ ■ ■ {fk^)wUm) 

with < ki < g and where Yl^iJi ~ 2r + 1. The further requirement that these terms 
are coefficients of the power n"'^^^^'' in Fr'^^ means that "^^iLi ki = 2g — 2r. Ehminating 
r from these two equations we may conclude that each term of Forcing^ is proportional to 
an expression of the form (|i5]) where ji + ■ ■ ■ + jm = '2{g — ki — k2 — ■ ■ ■ — km) + ^- This 
establishes statement (ii) of Proposition 13.11 



3.1.3. Laurent Expansions of Forcing g. In this section we set up and prove the main lemma 
needed to establish the pole expansions, Proposition 13.11 (iii). of the fundamental factors, 
(/fc),„(n)('S,l)- From ([27]) one has 

One can calculate the first couple derivatives directly: 

ifk)^ = (l-2A:)«;-2^z, + (z.-lK-i-2'^s4 
{fkU^ = -2k{l - 2k)w-^^''+^^Zk + l){iy - 4k)w''-^-^hz', + {u - 1)2^2.-3-2^^2^// 

The general form of the higher order tf-derivatives is given by 

n 



j=0 



Note that the coefficient of w^'^ {2k+n ^) gj j^^s two sources in (/fc)^(n-i). These two 
sources are 



and {v - lyP^- 



jp{n~l,k)^,(u-l)j^(2k+n-l)+l^j^(j) 



Given this, one can read off a recursion formula for the coefficients pj"''^^ (^u) that is given 
in the next lemma. 

Lemma 3.1. 

n 

(46) (A).(.)(s,l) = - ly Pj''^'\^)s^ zl^^ 

3=0 
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where 



(47) ''^V) = + {{^ - l)i -{2k -2 + n)}P,("-^''^)(i 

with 

(48) P(5"'^'^(z/) = (l-2A;)(-2A;)---(l-2A;-n + l) 

(49) Pk'^'^Xv) = 1 and 

(50) pj^'^'V) = for j > n and j < 0.. 



We next turn to studying the Laurent expansion of Zj. around zq = This analysis is 
subordinate to the inductive assumption that for < < g', 



2:0(^0-1)^3^-2(^0) 

^fe(^o) ^ 



{y-{y- \)z^f^-^ 



It fohows that 



(fc.O)/ N (fc,0)/ N „(fc,0) / N 



^k = Zo< — — + — — + 



{u-{u-l)zof'' {ly - {u - l)zof''+^ {u - {u - l)zof''-^ 



Moreover, differentiating this expression j times, with respect to s, one arrives at 



pk-i+j jk,j) 

(51) 4'^ = ci^rM E 



(^ - (Z. - l)^0)2fc+£+.- 
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To see this, note that 



~ ds 



^1/ ^0 ^0 / > 



^ (2fc + ^ + (j-l))a 



,-1 U-i)u , [(j-l)(i^-l)-2fc-£+l]af'^--^) 

+ V 2^ 



(z.-(^-l)zoP+^+^- 



= ^ .-+1 r's^^' [Q- - 1)^^ - (2fc + ^ + (j - 2)]af '^--^^ + K2fc + i + 2))af_'r^) ' 

where, in the fourth equation, the identity (j32]) . Zq = ^ = ^73^1:3^, was apphed. 
From this one deduces the following recursion for the succession of derived Laurent co- 
efficients 

Lemma 3.2. 

af''^ = [{j - l> -{2k + i + (j - 2))]af '^'-^^ + i^[2k + i + {j - 2)]af_'f 

with 

ar = 1 

af°) = for i>0 

(52) af '^'^ = for i < and i > 3k + j. 

Finally we consider the form of the expressions s^zj; which appear in (j46p . Recall, from 
(I5U1) . that 

(53) s = 
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SO one may deduce from ()3.1.3p that 

s^i^^ = zoizo-iy Yl 



{l-{iJ-{u-l)zo)y a} 



ZO 



l)j X] X] 



^ ' r=0 i=0 ^ ^ ' ^' 



ZO 



^ ^ r=0 m=r ^ ^ ' ^' 

where m = ^ + r. Inserting this last expansion into the expansion for (/fc)^(n) (s, 1) given 
by (|36]) . and applying the vanishing conditions ([52|) to extend the bounds of the inner 
summation, yields 

3fc+2j-i (-l)^-^' { I \ 



r 



ZO Yl 



m=0 



(54) = ^0 Yl 



In the second line we applied (j52p again to extend the upper bound of the innermost 
summation on the first line in order to be able to pull this summation to the outside. In 
the last line, the upper bound of the innermost summation was changed from j in the 
previous line to m. To justify this note that if m < j, then the terms with index greater 
than m in the original sum vanish anyway by (|52p : on the other hand, if m > j then 



terms greater than j in the new sum will vanish since ^ ^ ) = if r > j per the usual 
conventions for binomial coefficients. 



The proof of Proposition I3.1( iii) thus comes down to establishing the following vanishing 
lemma. 
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Lemma 3.3. 



n m / ■ \ 

(55) E(-^y~' i °— = 

j=0 r=0 ^ ^ 



for m = 0, 1, . . . n — 1. 



Proof. Fix A; > 0. The argument proceeds inductively. The base step is for n = 1 which 
imphes that m = 0. So the expression to consider is 



P(J^'''V) + '''V)(2fc - 1)) by Lemma 



[l-2k + l-{2k- 1)] ^''•"V) by dM 



For the induction step, we assume the lemma is true for n — 1; i.e., 



n— 1 m ^ \ 

j=0 r=0 ^ ^ 
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for m = 0, 1, ... n — 2, and consider the expressions on the left-hand sides of the putative 
equations in Lemma l3.3t 



j=0 r=0 



r I J 

r 



j=0 r=0 



j=0 



r=0 



n— 1 m 
\j p{n—l,k) ■ 



j=0 r=0 
n— 1 m 

j=0 



J. I m—r 



r=0 



+ (z^- l)E(-iy ^j" ^'''^^{-'^Yj ( ) a^~r by dSni) and shifting j in the middle sum 



j=0 r=0 

n— 1 m 



J, / m—r 



J + 1 

r 



{jV — (2/c + m — r + (j — l))}a, 



= -(2fe-2 + n)E(-l)^'^^""'''^E(-l) 

j=0 r=0 
n— 1 m 

+ iy{2k + m-r+{j- l)}a^:^]_i 

n— 1 m 

+ (--l)E(-l)''^^""'''^E(-irj'( i by Lemma EJ 



m—r 



j=0 



r=0 



= -(2fe-2 + n)E(-ir^^^""'''^E(-ir 

i=0 r=0 

n— 1 m 

+ (2fc-l+m)E(-ir^^^""'''^E(-l) 

j=0 r=0 
n— 1 m 



J, / "'m—r 



I^m-r ^""m-r-lj 



+ E(-i)''^^""'''^E(-i)''(^' 

j=0 r=0 
n-1 



{ 



j=0 r=0 ^ ^ 



n-1 



+ (z^-i)E(-i)^'^^""''T(-ir^' 

j=0 r=0 



J, / m—r' 
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We next rearrange the terms in the last expression and make use of Pascal's identity, 

r ) (r)~'~(r — 1 

to achieve some reductions. 

m 

,{n,k)X^f_,-^r I J \ (k,j) 
^m—r 



n—\ m / • \ 

= (m - (n - 1)) ;^(-l)^P^(-^''=) ^(-1)'- I o!^^ 

j=0 r=0 ^ ^ 

n— 1 in / . \ 

j=0 r=0 ^ ^ 

n— 1 m / • \ 



j=0 r=0 

n—1 rii ^ \ 

j=0 r=0 ^ ^ 

n—1 m ^ \ 

j=0 r=0 ^ ^ 

n—1 m ^ \ 

= (m - (n - 1)) 5^(-l)^P^(-^''=) ^(-1)^ aS:^] 

j=0 r=0 ^ ^ 

n—1 m / • \ 

j=0 r=l ^ 

n—1 m—l , 



^m—r 



r 

j=0 r=0 

n—1 m / ■ \ 

j=0 r=l ^ 

n—1 rn / • \ 

= (m - (n - 1)) ^(-1)^P,("-^'^) ^(-1)^ aS:^] 

j=0 r=0 ^ ^ 

n—1 ( m—l ^ . \ iii—^ 



j=0 r=0 
n—1 f m—l ^ . \ m— 2 



j=0 K r=0 ^ ^ r=0 
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The second line of this final expression vanishes by induction for all m < n — 1 since the 
upper bounds of the internal summations, m — 1 or m — 2, are < n — 2 in this range. 
Similarly, the first line vanishes, for m < n — 2, by induction. For m = n — 1 it also 
vanishes because the leading coefficient, m — (n — 1), vanishes for this value of m. This 
completes the induction and the proof of Lemma 13.31 □ 



4. Applications 



4.1. Integral and Recursion Formulas for Map Counts. In this section we will give 
several representations of the count for two-legged maps ()12p . The most direct one will 
involve the Laurent expansions of z^-" that were developed in Section [3.11 

We make use of the defining implicit relation, (jSOp . for zq given by 
(56) 1 = z{s) - az{sy 

where a = c^s and 




The j*'* coefficient of the Taylor expansion of Zg as a function of a near 0, Zg = J2j>o i^)'^'' ^ 
is naturally given by 



One can change variables from a near 0, in this integral, to z near 1 by differentiating the 
relation 

dz_ _ z" 
da 1 — vaz^'^ 

and then using the relation again to eliminate a in the differential 

da 1 — vaz'^^^ u — {y — l)z 

Iz ~ z^ ~ ' 

so that 

■(a),.., _ 1 / (•'-(i'-l)z)z'''-'z,{z) 



1 / .-p,,-.(.) 



2vri J^^i [z - -[v- l)zf9-'i 



32 N. M. ERCOLANI 

In the second line we have rewritten the integrand using the more expHcit form of Zg(z) 
given by Theorem 11.51 An alternative expression can be found by using instead the partial 
fractions expansion of Zg given by (j4ip . 



where X]i£o^ ^i'^^^^^ — since ^(-z) vanishes at z = 1 for 5 > 0. 



These integral formulas for the generating function coefficients may be re-expressed as 
recursion formulas. To accomplish this transformation we return to the a contour integral 
(j57p which may be recast in terms of higher derivatives as 

(58) = 



^^^^ ^ ^^L^ § 



where in the second line we applied the Laurent expansion (|51|) . 
Corollary 4.1. 

— *{two-legged gi-maps with j 2z^-valent vertices } = j\cPyC!'f\i') 



9 



4 af^\v) and 



£=0 

V) = [{3 - - (2<7 + ^ + {3 - 2))]a'-t^-'\v) + i.[25 + ^ + (i - 2)]4^jr^V) 
by Lemma 13.21 with 



,(0,0), 



1 



ai°'°)H = for ^>0 



,{9J)/ 



and 



for ^ < 



39-1 

£=0 
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The coefficients a\ ' (z^) appearing in this evaluation are thus recursively expressible in 
term of the 2>g coefficients, af'^\v), of the partial fractions expansion of Zg. These latter 
coefficients are rational functions of the , for k < g, over the rational numbers Q. 

4.2. Rationality of the Asymptotic Partition Function Coefficients. We will now 
make use of the following recursive integral formulae for the fundamental coefficients eg. 
(The derivation of these formulae is accomplished by a process similar to that described in 
section lA.21 for the coefficients Zg] for the details we refer the reader to |10j.) 

Theorem 4.2. [Tn] 



(60) 
egi-s) 



driverSg(2;o(s)) 



(2 -2(7) (1-2(7) ' 



1 



(s) _ ;l\ (29-2)/(i'-l) rzo{s) / ^^^u \ {2g-2)/iu-l) 



2 - 2(7 V CuZoisY 



z-l 



{drivers g{z))' dz 



+ 



1 



(1 - 2(7) V c^zq{sY 
when 5/1, where 



— — - I (driversg(2;))* (iz 

z 1 / 



(61) driverSg(2;o(-5)) 



^ (2£ + 2)!9?/;(2^+2) 

/ 00 ^ 

+ the n~^^ term of log I — 

Vm=o ^ 



W=l 



- Zm, [ S 1 



We denote by (driversg(z))* the derivative of diivevsg^z) with respect to z. K\ and K2 are 
constants of integration either determined by the requirement that Cg be a locally analytic 
function of s or by the evaluation of Cg for low values of v through its combinatorial 
characterization. When g = 1, 



ei(-s) 



1 



{s)-l\ f'ois) / c,z'^\ {u-{v- l)z) 



(62) 



(1/-I) \CuZq{s) 

Ms) (y-{y- \)z) 

z{,z-\) 
^log(i/- (i/- l)zo(s)), 



z-\ 
driversi(z)dz 



CyZ- 



-driversi(z)(i2; 



where we have chosen the principal branch of the logarithm. 
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For convenience we also record here the planar (g = 0) result, 
(63) eo{-s) = l\og{zo{s)) + ^^^^^{z,{s) - 1) (^zo{s) - 

We can now extend this characterization by proving Theorem 11.31 

For g >2, the coefficient eg{—s) = eg{zQ{s)) is a rational function of zq with poles only at 
zq = V l(y — 1) and which vanishes at least simply at zq = 1. 

Proof. We proceed inductively assuming the theorem to hold for all Cr with 2 < r < g 
and with ei and eo as given, respectively, by ([62]) and ([63|) . To start, one needs to examine 
the terms of driversg(z) for g > 2. Note that for 1 < i < g — 2, 



= (2 - 2{g - £)) {-w^-^) + w'^-^^^-i) {-w''~h) 

= (2 - 2ig - I)) ^i-2(^"^)e,_, {-w^-^s) + w^-'^^''^'-^el_, {-w'^-^s) 

CuZq 

where in the third line we have used the identity (j53p . The final expression for (|64p has 
three properties that will be useful to note: 

(65) It vanishes at zq = 1. 

(66) Its minimal pole order at zq = i//i/ — 1 is not less than that of eg-i{z). 

(67) It may also have poles at = but nowhere else. 

It is straightforward to see that these three properties are maintained under further differ- 
entiation with respect to w. 

We next separately check that these same properties hold for w-derivatives of eg-i when 
£ = g — 1 or g, which correspond respectively to ei and cq. For ei we have 

d r ^-l .1 _ W''-'^ Zq - 1 

dw L^'^ ""^J " 12 ' v-(v-\)zQ ■ czl ■ 

It is clear from this that the 2g^^ ^i;-derivative of e\ [—vf~^s) appearing in drivers^ has the 
three properties (|65l - I67p . Similarly one can see from (j63p that the {2g + 2)*^* it;-derivative 
of w^eo {—w^~^s) appearing in drivers^ has the three properties ([5S] -[S7 |) . Property ([55]) is 
vacuous in this case since the expression has no poles ?!X zq = v jv — \\^& take the minimal 
pole order in this case to be oo. 

Lastly we consider the last collection of terms in drivers^. 

1 Zrfni^s") \ 



( X] ^^™(^) 1 



log(zo(s)) + log 1 + 



oo 



m=l ^ ' 
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where Q3g-2r is a polynomial of degree 3^ — 2r. The form of the right hand side follows 
directly from Theorem 11.51 and shows that this part of drivers^ vanishes at zq = 1, has 
minimum pole order of 2^ at zq = -^jz^ and of order oo at zq = 0. We may thus conclude 
that 

Lemma 4.1. For g > 2, driversg(2;) vanishes at least simply at zq = 1, and its only 
singularities are poles restricted to occur at either z = or z = v l(y — 1). 

Also, driversg(z)' has only poles for singularities which are restricted to occur at either 
z = ^ OT z = V I {y — X) with minimal pole order at least 2 at both locations. 

We also observe, as was noted in [lOj, that the terms of (j6Up involving K\ and can only 
have the possibility to contribute if (2g — 2)/{v — 1) or {2g — l)/{u — l) are integers. But 
in that case it follows from (j53|) that they are terms which vanish simply at = 1 and 
have poles only at = 0. As a consequence of this observation and the lemma one may 
conclude from the representation ()60p that eg, for g' > 2, is locally rational at zq = ^j^. It 
thus suffices to consider the behavior of ([60]) near zq = 0, 1 in order to establish Theorem 



when s = (or equivalently when zq = 1) for all g. For the integral terms in (j60p . since 
by the preceding lemma driversg(2;)* is regular at z = 1, one may write the integrand in 
each of these terms as a Puiseux series around z = 1. The integral can then be carried 
out term by term on this series with constant of integration set to zero (note that the 
constants of integration are already captured in the terms containing the constants Ki 
and K2 ■ However, the prefactor in front of each of these integrals reduces the total 
result to a Laurent series at zq = 1 except for the possibility of one term proportional to 
{zq — l)'^log(2;o — !)• However, since by Theorem 11.11 Cj is analytic near zq = 1, such a 
term must either cancel with other such terms in the full expression for Cg, or it cannot 
exist. It follows that ()60p has a regular Laurent expansion at = 1 which in fact must be 
a Taylor expansion with vanishing constant term since eg = at zq = 1. 

We apply Theorem IB . 31 once more to conclude that eg is regular at zq = 0^ 



[L31 



Returning to Theorem 11.11 note that when t = 0, log 




) 



0; therefore, eg 



= 



N^9 N-^\og 



( 



ZN{t) 

Zn{0) 



) 



eo{t) 



eg-i{t) 





a 



36 



N. M. ERCOLANI 



uniformly as t — )• oo and hence Cg remains bounded as zq approaches 0. It fohows, in 
exactly the same manner as was done for zq near 1, that all singular parts of (j60p at = 
must cancel each other. □ 

Corollary 4.3. The integration constants Ki and K2 can be non-zero only if {2g—2) / (i^—l) 
or {2g — — 1) are integers and in addition the combination of the other terms in (j60p 
has a pole at zq = 0. 

4.3. The Double Scaling Limit for u = 2 and Painleve I. The topic of the so- 
called double- scaling limit of the random matrix partition function ([3]) has been extensively 
discussed in both the physics and the mathematical physics literature related to two- 
dimensional quantum gravity. For a concise history of these developments and relevant 
references we refer the reader to [14]. In this literature attention was focused on the 
special case of the weight ^(A) = ^A^ + tX^ which corresponds to the case of = 2 in 
our more general treatment. In this case the double-scaling limit refers to the regime in 
parameter space where t — )• tc{= ^) and N ^ 00 simultaneously so that 



remains fixed. The authors in [31 El IS] arrived at this prescription of the limit through a 
formal scaling argument based on the Freud equation (jSlh . Basically they posited a general 
scaling form, in n, N and t, for the recurrence coefficients 6^ ^(t) and determined that ()68p 
was the unique scaling that would (formally) asymptotically balance (jSip . Realizing that 
this relation, which they called the string equation, could be taken to be a discrete form 
of the first Painleve equation 



they took ([69]) to represent the exact non-perturbative string theory for two-dimensional 
euclidean quantum gravity. 

The corresponding mathematical conjecture was that the large N expansion of the recur- 
rence coefficients, under the scaling (j68p . has a limiting form y{^), at leading order in 
A^, that solves (j69p which hereafter will be referred to as PI. The proposition was first 
placed on a rigorous footing, based on Riemann-Hilbert analysis, in the papers fl2l I13j. A 
complete and self-contained proof was recently detailed in |7|. 

The Freud equation and its relation to PI play no role in our description of the double 
scaling limit for general v which we take up in the next subsection. The reason that we 
have dwelled on them for the case = 2 in this subsection, apart from their historical 
relevance for our topic, is that this special case provides a concrete point of comparison 
(see Section 14. 5p for our more general results and a jumping off point for possible future 
investigations that might involve a broader connection to Painleve transcendents. For these 
reasons we will briefly review, here, the essential properties of the relevant PI solutions. 
For a systematic development and further details we refer the reader to [19j . 





(69) 



y" = + e, 
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There is a family of solutions to ()69p characterized by having the common asymptotic 
expansion 




(70) y{0 ~ Vlf h + 



-oo 



for particular coefficients Og. These expansions are valid in angular sectors, typically of 
width ^ in the complex ^ plane. They differ from one another by jumps across sectors 
(Stokes lines) which are exponentially small. When needed, for specificity, we will take 
y(^) to be that solution with the above asymptotic expansion valid in the full sector ^ < 
arg^ < as |^| — )• oo. The ag satisfy the following quadratic recursion formula 



^7^^ 25g^ - 1 



m= 

ao = 1. 



4.4. The Double Scaling Limit for General v. In this section we will describe a 
double scaling limit for the free energy of the unitary ensembles ([T|) in terms of the rational 
representation of the Zg . This will be done in the setting of general values of u which should 
enable one to see the universal character of these results. 

We return to the asymptotic expansion for the orthogonal polynomial recursion coefficients 
(Theorem 11.41 or Theorem IB.3P and substitute the polar expansions found in Theorem 
[131 



(72)6^,^ = zo + £%iV-'^ 




oo 

2<7 



9=1 



(l/ - (j/ - l)zo)29 (i. _ (z, _ l)^o)2g+l _ (j, _ l)2Q)5g-l 



We first proceed formally to determine what the double scaling balance should be; i.e., we 
choose 6 in 

l)zo ~ 

such that the highest order terms of (I73p at the pole have a common factor involving N 
that is independent of g. In other words 6 should be chosen so that \_{5g — 1)6\ = —2g. 
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Therefore, 5 = — |. Combining this with one is led to take the double scaling ansatz 
for general to be 

(74) 



where, by (j33p . Sc = ^'^ ^ — and ^['^'^ is a constant depending only on u. We will in general 
notationally suppress the explicit dependence of 71 on i' since this should be clear from 
context. Also, by 



v-l 



S-Sc){l + 0(S-Sc)} 



-71 



Substituting the ansatz ([7^ into ([73]) , one may formally conclude that 



{u-iy 



-71 



e Ar-2/5 



9=1 



+ ••• + 



71? 



(5g-l)/2 



-2ff 



V (^-1) 



Ar-2/5 



(75) 



+C)(Ar-4/5 

-2c;,7ii^ 



1 '^^ u-l^_ 



(53-l)/2 



\l/2 



(zy - 1)^+2 



00 



9=1 



-20^711^" 



1 \ -59/2 



Ar-2/5 



We next study what can be said about the leading order, for large N, term in the formal 
expansion (|75p . The coefficients of the series at this leading order are essentially given by 
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the a^g_i{i') which are, sequentially, the highest order coefficients of the pole at — 1) 
for the sequence of Zg. We derive a recursion formula for these coefficients. 



Proposition 4.4. 

(76) oSd-iI") = 

4'V) = 



1/3 (255^ 



(m) 
'''3m- 1 



f \ (g~m+l) I N 
^'')«3(g-m+l)-l(^) 



m=l 



(77) 



Proof. First recall from the proof of Theorem 11.51 given in Section [3] that the coefficient, 
that a3^^_]^(z^), has (inductively) two sources. The first of these comes from forcing the 
integral equation for Zg (see Theorem 1 2 . 41 ( iii) ) with terms of the form 

9-1 

(the J = 1 case) . The other comes from forcing with terms of the form 

( f ) 

6J \Jg~^Jwww 
J w=l 

(the J = 3 case; note that here we have used Proposition 12.21 to calculate that d^^-^^ = 

Inserting the partial fractions expansion of these terms as the forcing in the integral equa- 
tion for Zg and then collecting the coefficients of the terms proportional to l/(z^ — (i/ — 

1)zq)^^~^ one derives an equation involving a mixture of coefficients of the form af''^\v). 
However, by repeated application of the linear recursions from Lemma 13.21 supplemented 
by (j48p - (jSOp . one arrives at the following expression involving only the maximal pole 



coefficients of undifferentiated Zr, 



9-1 



g-m 
3{g-m)- 



m=l 



+ 



6(i/ - 1) 3(3- 

Making some obvious cancellations and shifting from g to g + 1 one thus deduces (j76p . The 
initial condition ()77p follows directly from examination of (I13p . □ 

Using (j77P one can rewrite ([7U|) as 



«%+i)-i(^) 



It is immediate from this that a 



6 

(9) 
39-1 



9-1 

\ ^ (m) 



(g-m+l) 



"'3m-ll^J'^3(c,-m+l) 



m=2 

'z^) > and hence. 
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Corollary 4.5. The pole order of Zg at zq = v/{y — 1) is exactly equal to 5(7 — 1 for all 
It also follows from this form of the recursion that 

and therefore the series in the leading order term of ()75p is divergent as ^ — t- — oo. One 
expect this series to be an asymptotic and their are various approaches one could take 
to determine to what functions (as v varies) it might be asymptotic (generalized Borel 
summation, seeking an ode with a solution whose asymptotic coefficients satisfy Proposition 
14.41 etc.). However, as we will suggest in the next section, the form of our recursion suggests 
another approach both to characterizing the leading order series and to establishing the 
validity of the double scaling expansion ([75ll . 



4.5. Universality and the Painleve I Hierarchy. We will now restrict our study of 
to the special case u = 2 where it specializes to the following, still formal, expansion: 



(78) = - (-19271^) 



1/2 



4?-i (-1927iO"^^^' > N-^'^ + 0{N-^'^) 



We quote a recent result of Duits and Kuijlaars which will enable us to give a precise 
characterization of (1781). 



Theorem 4.6. [7] There are constants 71 -see and 72 such that 

(79) bl^j, - 2 = -72 (y„(e) + yMO) N-^"' + ©(iV"'/') as iV ^ 00, 

and where ya and are two members of the family of PI solutions having the common 



asymptotics (70). The constants 71 and '^2 o,re independent of the choice of a and p. This 
expansion holds uniformly for ^ in compact subsets of R not containing any of the poles of 
Ha and yp and can in fact be extended to a full asymptotic expansion in powers N~'^^^ . 

The derivation of the above result depends fundamentally on an extension of the defini- 
tion of orthogonal polynomials with exponential weights to a more general class of non- 
Hermitean orthogonal polynomials which corresponds, in an appropriate sense, to taking 
t in the weight V to lie in the interval ("fg-,0). There is a corresponding extended notion 
of the equilibrium measure mentoned in Section 12.51 for these negative values of t. In [7] 
this measure is defined through a variational problem along a deformed contour in the 
complex A plane. These orthogonal polynomials are characterized by a Riemann-Hilbert 
(R-H) problem and their asymptotics may be studied via the Deift-Zhou method of nonlin- 
ear steepest descent [6]. In particular, this involves the construction of local parametrices 
in terms of the Riemann-Hilbert problem for the PI equation. (The parameters a and /3 
in the above Theorem extend through this process from the contour deformations to the 
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R-H problem for the orthogonal polynomials to the R-H problem for PI to the family of PI 
solutions. A key point for relating this analysis to our analysis in the case of (j78p is that 
for ^ < t < there is an NQ{t) such that for all N > N(j{t) the recurrence coefficients 
b% ]\i{t) for the non-Hermitean polynomials exist ([7j, Theorem 1.1). Moreover, b% ]\i{t) 
has a full asymptotic expansion, independent of the choice of a and /3, in inverse powers 
of N. Indeed the construction of this expansion is modeled on that of [9]. It follows that 
b%,Ni*) = ''N,Nit) for ^ < t < 0. 

We can use these observations to directly relate the highest order polar expansion coef- 
ficients, «3pLi of Zg for z/ = 2 that appear in ([78|) to the coefficients in the asymptotic 
expansion (|70|) of the PI solution. First we need to determine the pinning constants 71 and 
72. Matching the coefficients of N^"^/^ between ([75|) . ([7^ and ([70]) requires 

-(-192716'/' I l-2f;ag_i(-1927iO-'^/4 = -I2{ya{0+ym) 



9=1 



(80) = -272V4|l + E«5(-0 



-55/2 

9=1 



Comparison of the first two terms gives 



4 9 

19271 = -72 



(19271)5/2 



6 



ai = -^by^n 



^2 



\ by m- 



From this one immediately sees that 

(81) 72 = 2^/532/5 

n-9/5o-6/5 1 

(82) TI = = 

These values agree with the pinning constants stated in [7] (Theorem 1.2) modulo the 
factor of 1/4 in ()82p . However, in that reference the form of the weight is taken to be 
V = \^ /2 + tA^/4 which differs from the form of the weight used here effectively by scaling 
the time t by 1/4; hence, the first pinning constant defined in (j68p should differ in fact 
from that appearing in [7] by a factor of 1/4. 



With these constants determined one may now use (j80p to express the highest order pole 
coefficients (for u = 2) m. terms of the PI asymptotic coefficients. 
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a. 



(g) 



25S-1 (2/3)f/2 for g > 1. 



Substitiuting (f83]l into the nonlinear recusrsion (fTTI) one deduces the following quadratic 
recursion between the coefficients of the highest order poles of the Zg in the case of = 



This is in complete agreement with Proposition 14.41 for the case u = 2 and so we can finally 
state 

Proposition 4.8. The leading order series in (I78p coincides with the asymptotic expansion 
of the PI solution specified in Theorem l4.6l Moreover, based on the validity of the expansion 
(f79l) . the expansion (|78l) is also valid for large N. 

In addition, from the derivation of the expansion displayed in ()75|] one deduces an improve- 
ment of the result stated in Theorem 14.61 

Corollary 4.9. The double-scaling limit of bf^ has a full asymptotic expansion in even 
powers of N~^^^. In particular, the next order correction in (I79p . after leading order, is 



The form of the expansion (1751) suggests that Proposition 14.81 for the special case of = 2, 
should extend to have a universal character for general i'. By universal here we mean 
that the various constants and expressions we have been discussing when v = 2 should 
be replaced in the general case by fixed rational functions of the parameter i/. Although 
we do not carry out the details here, this extension should follows straightforwardly by 
generalizing the Riemann-Hilbert problem for PI to the R-H problem for the PI hierarchy. 
(For a description of this hierarchy we refer the reader to [26|.) Specifically one needs to 
replace the Hamiltonian for PI appearing in the RH problem by the Hamiltonian for the 
higher order equations in the PI hierarchy. 

Another intriguing problem is to study the fine structure of the higher order terms in 
the double scaling expansion (|79p ; i.e., the coefficients of N'"^^/^ for integer h > 1. This 
should lead to extensions of Proposition 14.41 that enable one to recursively determine the 
fundamental coefficients af\v) for lower values of I.. 

4.6. Relation to Other Enumerative Results. The subject of map enumeration has 
been much studied going back at least to the early work of Tutte [28], in the '60s, who 
introduced the notion of a rooted map. A map is said to be rooted if a vertex of the map 
together with an edge adjacent to it and a side of that edge are distinguished. Since then, 
much work been done on counting various classes of rooted planar {g = 0) maps. The cases 



2. 




m=l 
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of higher genus maps have been more challenging and therefore the main emphasis has been 
placed on trying to determine the asymptotic behavior of various types of enumerations 
for large values of the discrete parameters.. A particular example of this is given by the 
problem of counting the number, Mn^g, of rooted maps on a genus g Riemann surface with 
exactly n edges. Bender and Canfield [Tj showed 

(84) M„,g ~ tgn^^3-i)f2 ^ n^oo. 

Many other map classes have similar asymptotics [15j with the same sequence of constants 
tg. In ^ a recursion is given for these constants the first few of which are: 

2 1 7 



24 ^/4320 



IT 



However, there have been some recent improvements in the recursion formula for the tg [2] 
based on connections to the enumeration of branched coverings of Riemann surfaces [17] . 
More recently [16] it has been observed that this new recursion is closely related to the 
recursion for PI asymptotic coefficients (|7ip . The upshot is that 

(85) tg = -^-^^fO 

= ^ ■ f(5k)«3?-i by m, 

which reveals an interesting albeit mysterious link between two prima facie quite different 
classes of map enumeration problems. 



Appendix A. Recurrence Coefficients and their Continuum Limits 



This appendix provides a reasonably self-contained presentation of the relevant results 
(and outline of their proofs) from [Qj [10] . It also gives the proof of Proposition 12. 2t a new 
characterization of the coefficients in the continuum Toda equations. 



A.l. Full Asymptotics of the Recurrence Coefficients. Recall from Section 11.4.41 
that we have defined 7r„,Ar(A,t) to be the n*^ monic orthogonal polynomial with respect 
to the exponential weight exp(— A^y(A)) and that these polynomials satisfy a three-term 
recurrence relation of the form 

(86) 7r„+i_Ar(A) = A7r„,Ar(A) - 6^_jv(0^n-i,A'(^)- 

There is a basic relation between these recurrence coefficients and the tan functions defined 
in ([6|) given by: 
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(87) bl,{d) = i^log[r2, (01,0)]^^^^ where 



Tl.iOue) = zi")(0i,^)/Z;")(O,O) and 
(89) 4"^(ti,t) = |---|exp|-Arf;QA2 + af + tiA,Uv(A)d"A, 



where V( A) = JJ|Aj-A/p 

j<i 

The relation ()87p is called a Hirota formula [11] and and stems from the integrable systems 
theory associated to the Toda Lattice (see Section IA.2P . The connection to orthogonal 
polynomials comes from the representation of the partition functions ^ as Hankel deter- 
minants of moments for exponential weights [9l \TU[ [27] . The definitions on the next two 
lines are consistent extensions of our earlier definitions. Indeed when Oi,ti = 0, (jSSp and 
([89]) are equivalent to ([6]) and ([3]) respectively. 

The Hirota formula may be extended to more generally scaled recurrence coefficients by 
making use of the scalings introduced in Section I2.1.1[ Indeed, the tau functions are 
convariant with respect to these scalings so that 

(90) Tn,i{0i,9) = Tn.N{ti,t) = T„,„(— Si , —s) where wc define 

(91) ti = 291 /Vn consistent with and 

(92) si = —ti/y/x consistent with (fT8]l . 

The final transformation to is realized by a variable change A = ^/x\ of the eigenvalues 
in (IMl). It then follows from (fel) that 



1 (P 

= [rlni-si, -s)] ^^^0 ; moreover, 

N 1 — s 1 (f 

(93) _62^^(t) = _62_^(__^) = __log[r2„(-si,-s)]^^^g. 

Theorem 1 1 . 1 1 mav now be applied to the right hand side of ([93|) to establish the expansion 



A. 2. The Continuum Limit of the Toda Lattice Equations for the Recurrence 
Coefficients. As the external weight, VJ/(A;t) changes with t, the corresponding orthogo- 
nal polynomials, vr„^7v(A, t), also evolve. How they evolve is governed by how the recurrence 
coefficients, which also depend on t, change. It is well-known that the t dependence of the 
recurrence coefficients is governed [4j by a system of nonlinear differential equations known 
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as the Toda Lattice Equations. There is a different system for each z^. Though different, 
these different nonhnear systems have a common general form. The different associated 
flows on the recurrence coefficients that these different systems induce commute with one 
another and, in fact, this commutativity is related to the complete integrability of the Toda 
Lattices. This feature does not play a direct role in what we present in this paper but it 
is undoubtedly related to the "universality in z^" that has already been mentioned several 
contexts. 

In our setting, the Toda Lattice equations take the following form, for each v. 



(94) 



2 cie 



E 



.m=l 



m=l 



n+ern(w) 



where the dependent variable 5^ stands for 6^ i{9) as defined in ([87|) . (A detailed deriva- 
tion of the form (j94p of the Toda equations from their more standard presentation in the 
integrable systems literature is given in Section 4.1 of [ID].) The summation in (j94p is 
taken over the set of walks, denoted {w}, of length 2z^ on Z that start at +1 and end at 
— 1. If one visualizes the Z lattice along which the walk takes place as a vertical axis and 
the steps of the walk as discrete equally spaced points along a horizontal axis which serves 
to order those steps, then the walk is graphed as a zig-zag path comprised of line segments 
of slope +45° (upturn) or —45° (downturn). Such a walk is completely determined by 
specifying where its downturns (of which there are exactly v + 1) occur. If the discrete 
variable i denotes locations on the vertical axis then, on a given walk, i can only range 
over the interval [— i^, v]. Finally we let Imiw) denote the vertical axis location on the walk 
w after its m''^ downturn; then im{w) + 1 denotes the location before this downturn. 

IQl). that these 



In order to study the continuum limit of (j94p we note, from (j93p and 
equations contain terms with prima facie different asymptotic gauges: 



(95) hi = hl , 



?)2 - /,2 



1 1 

n [ zo{s) + -^zi{s) + -^Z2{s) + 



2(n + . 



(n + £) zq{s) + 



r^i(s) + 



-Z2{s) + 



(96) 



n [wz^{sw'-^) + • • • + ^y-''Hg{sur-^) + • • • ^ 



where w = 1 + - as defined in ()19p . The introduction of the latttice scaling variable 
w in (|96p allows us to analyze all the terms appearing in (j94p with respect to the same 
asymptotic gauge. Since I G \—v,v\ and v is fixed, one may assume that i « n. Then 
one can study the limit as n — ?■ oo of the equations (I94p . with the substitutions (I96p . as one 
would study the continuum limit of a numerical scheme or a molecular chain with analytic 



46 



N. M. ERCOLANI 



potential. To carry this out, we introduce a more compact notation for (j96p as 



(97) 



(98) 



f{s,w) 

fg{s,w) 
}? 



2 d9 



1 1 

fo{s,w) + ;:;2-^i(^'^) ^ T^fai^^'^^ ' where 



w^'"^^ Zg{w'^~^ s) so that 



n 



nf I s, 1 H and in particular 

f{s, 1) = Zo{s) + -^Zl{s) + ^Z2(s) + 



and 



1 1 ci5^ V ft N 



UI = 1 



Finally, substituting ([971) and dMl) into ([MI) and Taylor expanding in ^ (and ^^^) one 
arrives at the basic continuum equations 

E"/"*' n 



w=l 



{w} 

(99) 



V. m=l 

n 

m=l 



f \ n ) 2/ I n 



+ ••• + 



h\f 



^ / n + 2/ V " 



_|_ . . . _|_ ( ^ I _|_ 



hlf V n 



from which Theorems 12.11 and 12.41 can be deduced. 



n 



+ 



A. 3. Coefficient Formulae for the Continuum Toda Equations. We now provide a 
proof of Proposition 12. 2[ Recall from Theorem 12.11 that . using multi-index notation. 



9 



E 

|A| =25 + 1 
£(A) <u + l 



n,-^.(A) 



-1 1 g'^l/ 

A! dw^ 



These coefficients have the character of correlation functions for certain tied random walks. 
A description based on this perspective is presented in jlOj . In the following derivation, 
that description will be used but not explained in full detail. 



Proof of Proposition 12.21 The coefficients d^'^' are directly related to the expression of 
the Toda Lattice equations for the 6| , in terms of tied walks on a one dimensional lattice as 
described in [10]. The walks in question can be reduced to considering walks of length 2i^ 
on Z that start at +1 and end at —1. If one visualizes the Z lattice along which the walk 
takes place as a vertical axis and the time steps as discrete equally spaced points along a 
horizontal axis, then the walk is graphed as a zig-zag path comprised of line segments of 
slope +45° (upturn) or —45° (downturn). We label the locations on the vertical axis by 
i and those on the horizontal axis by i. For the walks under consideration, i ranges over 
the interval [—i^, i'] and i ranges from 1 to 2z^. Such a walk is completely determined by 
specifying when its downturns (of which there are exactly ly + 1) occur. We will let ij, for 
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j € {1, . . .u + 1}, denote the step at which the j*^ downturn takes place. We also let ij 
denote the location of the path after the j*^ downturn (then £j + 1 denotes the location of 
the path before the j^^ downturn). 

Given this set-up, and based on the derivation given in }10j, the coefficients can be expressed 
as 

dx''^ = E [mx{h + l,...,iu+i + l)-mx{eu---,i.+i)] 
walks 

(100) 

= ^ [mA (ii,i2 - 2, . . . , v+i - 2i/) - mA (ii - l,i2 - 3, . . . , v+i - - 1)] . 

ii<---<ii,+i 

We now re-express the sequence of downturn steps as a restricted partition as follows 

fJ, = {fXl, . . . , fXj, . . . , flu+l) = {iu+l, ■ ■ ■ ,'i'u+l-{j~l), ■ ■ ■ :h)- 

We also define the partition 

r] = (2i^, 2i/ - 2, . . . , 2, 0) 
and for notational convenience we define the difference 

fJ'-V = (wi - 2z^, V - 2z^ + 2, . . . ,^2 - 2,ii). 
Since mx is a symmetric polynomial we may rewrite (jlOOp as 

dx''^^ = Yl [mx{fi-v)-mx{f^-r]-{l,...,l))] 



'^*='^(i/+l,i/,...,2,l) 



where 



^('+M,.'.',2';it = {/u e 71 : (z. + 1, z., . . . , 2, 1) C ^ C {2v, 2u - I, . . . ,u)} . 



Note further that 

(_(j._l),_i.,...,0,l) C /i-r? C (0,1,. ..,1.-1,1.) 

(_i,,_i. + l,...,_l,0) C ^-7?-(l,l,...,l) C (-l,0,...,i/-2,z.-l). 
It follows that /i — 77 — (1, 1, . . . , 1) ranges over the same set of sequences as ii — rj, but with 

all signs reversed. Hence, the sum in the last expression for S^'^^ may be rearranged to 
read 

dx''^^ = Yl [mxil^-v) -mx{-{^^-r]))] 

„(2>.,2^-l,...,^) 
(■^+l,'^,-..,2,l) 



TTiA (/X - 7/) - (-l)l-^lmA (/i - r/) 



„(2^,2,.-l,...,^) 
(■^+1,1',. ..,2,1) 



Y 2mx{n-v), 



^{2,.,2,.-l,...,.) 
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where the last hne follows because |A| = 2g( + 1 which is odd. □ 

Remark 9. The formula for d^y^'^^ in Proposition 12. 21 differs significantly from that given in 
|10j . Both versions have their respective merits. The one introduced here is, for example, 
useful for growth estimates of these coefficients with respect to their parameters. 

Appendix B. Large-Time Behavior of the One-Point Correlation Function 

In this section we are going to discuss the large N behavior of the one-point correlation 
function for eigenvalues oi N x N random hermetian matrices in the vicinity of exponential 
weights Vt{\) = Vj/(A;t) for which t is very large. We refer to this as large-time behavior 
because t has the interpretation of a dynamic variable of the 2v^^ Toda flow. This discussion 
applies to the determination of the large-time behavior of the coefficients Zg and eg. The 
analysis of the large N behavior of the one-point correlation function that was carried out 
in [9j was done explicitly only for values of t close to zero. However, as was pointed out in 
that paper, the range of validity of the large N expansion derived there was by no means 
restricted to a small neighborhood of t = 0. Indeed, the expansion can be constructed with 
uniform validity in an open complex neighborhood of any t > as long as the equilibrium 
measure /ivi for the weight Vt at that value of t, satisfies the following conditions: 

(a) The support of ^Vt is a single interval of the form [— /3,/3]; 

(b) ^Vt satisfies the variational equations of (f37|) : 

(c) ^Vt vanishes like a square root at both endpoints of its support. 

In Section IB. II we will show that these conditions are indeed satisfied for all t > by 
explicitly writing down the equilibrium measure for all these values of t. It then follows, 
from analytic continuation and uniqueness of power series expansions for analytic functions, 
that the coefficients of the large expansions of logarithms of tau functions ([6]) and their 
derivatives, at these values of t, coincide with those constructed near t = (such as eg and 



However, for our applications we require more. We need to know the limiting values of 
these coefficients as t — )• oo. To accomplish this one should actually build the asymptotic 
expansions of correlation functions in a neighborhood of t = oo. In section IB. 21 we indicate 
how the constructions of [9] may be adapted to this large-time regime. 

B.l. The Equilibrium Measure for t > 0. We present here an explicit expression for 
the equilibrium measure (I35p that is uniformly valid for all t > 0. To accomplish this, 
we rescale the domain variable A so that the support of the measure remains fixed as t 
varies: 




(101) 
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Proposition B.l. The equilibrium measure for t > is given, under the variable dilation 
(jlOip . by the following expression which is a convex linear combination, for zq S [0, 1], of 
two probability measures on [—1,1]: 



dfMVt in) 



IT 



Zo + {I- Zq) 



(2r/) 



2u-2 



'2v-l 
.u-l 



(102) 



2i-l 



V(?7 + l)(l-r/) dri. 



This expression depends on t on/?/ through the parameter zq. When zq = 1 this corresponds 
to the weight for t = and the measure has the density of the Wigner semicircle law; when 
Zo = this corresponds to t = oo which we think of as the limiting asymptotic measure 
for large t. This limiting density also gives a probability measure supported on [—1, 1]. It 
is in fact the density of the equilibrium measure for the pure monomial weight 



(103) 



1?^ 

VM = —rfr 



Proof. Pinning down an explicit expression for ()35p amounts to determining the coeffi- 
cients of the polynomial h{\). We already know /3, from (|36p . explicitly in terms of zq. 
The essential ingredients for calculating 



(104) 



h{\) = i + ^/i^A^J' 



were developed in section 3 of |10j . We recall those elements here. 
Define the sequence {vjlj^g 



j=0 



whose Taylor coefficients can be computed to be 



(105) 



Vi = -7 



1 f2i - 1\ /3 



4n ^ - 1 J i + l 



for i > with vq defined to be /3^/2. The coefficients of h are directly expressed in terms 
of the {vj} as 
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(106) h, = 4u{u-j)t^^^ 



/32 

ZQ Vu-l-j 
/32 



K-1 



2v - 2j - 3\ 

) jl-zo) 

2.-1^ 

1 (l-^o) 
2v - 1 



for j < u — 1] 



Equation (j30p was used to eliminate t{= —s) and was replaced by its definition from Q. 
The exact expression (jl02p now follows by substituting (jl04p with the coefficients ()106p 
into ([35]) . using ([36|) to re-express /? in terms of changing variables according to (jlOip 
and, finally, changing the index of summation from j to v — j. 

It is an exercise to check that ^Voo is a probability measure by directly integrating the 
density (jl02p . for zq = 0, over [—1, 1]. (It is manifest that this density is positive on (—1, 1).) 
The equilibrium measure for monomial weights such as (jlOSp is explicitly calculated on page 
183 of [4J. One may check, by direct comparison, that (jl02p . with zq = 0, is indeed the 
equilibrium measure for (jl03p . The next section will give another perspective on why this 
is the case. □ 



B.2. The Riemann-Hilbert Problem at Infinity. The prior results (Theorems II. 1|, 
11.41 12.ip cited in Sections [1] and [2] , and on which the results of this paper fundamentally 
depend, all derive from a detailed asymptotic analysis of the spectral density for eigenvalues 
of random Hermetian matrices that was carried out in Explicitly this spectral measure 
is defined as the expectation 

{m)pi\t,\) = -^E^,, ij^*^^ • ^ (-OO'^)}) ^-^-t- tlie probability density 



1 ^ 1 



This has a remarkable expression as a one-point correlation function in terms of the Wron- 
skian determinant 

(108) pS^)(a) = kn{\a) = [yiiimiW - Y^,{\)Yi,{\)\ 

of the following Riemann-Hilbert (RH) problem for the 2x2 matrix Y{X): 
• Y analytic in 
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Y ={I+ A^'^^^ 

Y has Holder continuous boundary values Y±, from above and below respectively, 
along A G M, 

y+ = y_(i--7*'''),for AG 



where CT3 is the diagonal Pauli matrix referred to by this standard notation. This problem is 
directly related to a RH problem for orthonormal polynomials with weight Vj. We will not 
need to go into that here, but refer the reader to [9] for more details and references. 

The relevance of the one-point function to the problems discussed in this paper can be seen 
from the fact that 

(109) Z7v(t) = ^iv(0)exp j-iV^^ J p[^\t, X) X^^dXclt^. 

This falls within the more general class of problems which study the large N limits of 
integrals of the form 

(110) [ p(^)(t,A)/(A)dA 



where /(A) is a general C°° with no worse than polynomial growth at infinity. 

We now want to study how the RH problem for Y behaves in the limit as t — )■ 00. Under 
the change of variables (jlOip the problem for Y transforms to an equivalent problem 



• Z analytic in 

. Z = (/+ 0(i))r/^'^3, 

• Z has Holder continuous boundary values Z± along ?] G M, 

• Z+ = Z_ 1 e-^ytM ^ , for A G M, 
where 

(111) = 27?2 + (1 - zo) —ifr 

C-v 

Equation ([30|) was used to eliminate t(= —s). We see that as t — t- 00, — >• and the weight 
limits to '^ff'^^ consistent with (jl03p . This establishes that the equilibrium measures (I35p . 
for the variational problems (j37p . converge as t — )• 00 to the equilibrium measure for the 
variational problem with monomial weight (jlOSp that is associated to the RH problem for 
Z at zq = 0. Given this, all of the analysis of [9j to derive the full asymptotic expansion 
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for p\ uniformly valid for positive t near 0, carries over mutatis mutandis to give the 
asymptotic expansion of this one-point function near t = oo: 

Theorem B.2. There is a Zq £ [0, 1] such that for all zq G [0, Zq), the following asymptotic 
expansion is valid uniformly in zq: 

/oo 
fiX)p[''\tizo),X)dX = fo + N-^fi + N~^f2 + --- 
-oo 

provided that the function /(A) is C°° and grows no faster than a polynomial as |A| — ?■ oo. 
The coefficients depend analytically on t for ZQ{—t) G [0,Zq), and the asymptotic expansion 
may be differentiated term by term. 

This construction is reminiscent of the small amplitude limit for action-angle variables 
in Hamiltonian mechanics. It is, in effect, a one-point compactification at infinity of the 
family of variational problems ([37|) . their respective unique solutions, and associated RH 
problems all parametrized by t € [0, oo). This parametrization is now naturally referenced 
to the compact interval < zq < 1. 

Finally we apply this result to extend the domain of uniform validity in t for Theorems 11.11 
andO 

Theorem B.3. There is a constant A > such that for (complex) t with > 0, 

< A one has a uniformly valid asymptotic expansion 

(112) iogrf^^^{t) = Nho{t) + ei{t) + ^e2{t) + ■ ■ ■ 

as N ^ oo. Also, the recurrence coefficients for the monic orthogonal polynomials with 
weight exp(— A^y(A)) have a full asymptotic expansion, uniformly valid for (complex) t 
with 3f?(t) > 0, |9(t)| < A , of the form 

(113) b%^^it) = zoi-t) + ^zii-t) + ^Z2{-t) + --- 

as N ^ CO. The meaning of these expansions is: if you keep terms up to order N^'^^ , the 
error term is bounded by C N""^^""^ , where the constant C is independent oft in the domain 
{(Kt > 0; —A < 3=* < A}. Moreover, in this domain, for each I, the functions ei{t) and 
Zi{—t) are analytic functions oft and the asymptotic expansion of derivatives of log {Zj\f{t)) 
and b'ff j^(t) may be calculated via term-by-term differentiation of the series. 

Proof. For any t £ [0, oo) the methods of [9] enable one to construct asymptotic expansions 
of the form (I112p and (I113p which are uniformly valid in a complex neighborhood of that t 
intersected with the half plane JRt > 0. By Theorem IB. 21 we also have such a neighborhood 
around t = oo. Consider an open covering of [0, oo] by such complex neighborhoods which 
corresponds to an open cover of zq G [0, 1]. By compactness of this latter interval there is 
a finite sub-cover which in turn corresponds to a finite open cover of t G [0, oo) for which 
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the genus expansion ([7]) and the recurrence coefficient expansion (jlOp is uniformly vahd. 
Since the cover is finite, a width A in Theorem IB.3I can be determined. □ 

The uniform vahdity of these expansions in a semi-infinite strip of the t-plane is important 
for the results of this paper. See, in particular, section [3.1.1l and the end of section HT2l 
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